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We present a nonlinear post-Friedmann framework for structure formation, generalizing to cos¬ 
mology the weak-field (post-Minkowskian) approximation, unifying the treatment of small and large 
scales. We consider a universe filled with a pressureless fluid and a cosmological constant A, the 
theory of gravity is Einstein’s general relativity and the background is the standard flat ACDM 
cosmological model. We expand the metric and the energy-momentum tensor in powers of 1/c, 
keeping the matter density and peculiar velocity as exact fundamental variables. We assume the 
Poisson gauge, including scalar and tensor modes up to 1/c 4 order and vector modes up to 1/c 5 
terms. Through a redefinition of the scalar potentials as a resummation of the metric contributions 
at different orders, we obtain a complete set of nonlinear equations, providing a unified framework 
to study structure formation from small to superhorizon scales, from the nonlinear Newtonian to 
the linear relativistic regime. We explicitly show the validity of our scheme in the two limits: at 
leading order we recover the fully nonlinear equations of Newtonian cosmology; when linearized, our 
equations become those for scalar and vector modes of first-order relativistic perturbation theory 
in the Poisson gauge. Tensor modes are nondynamical at the 1/c 4 order we consider (gravitational 
waves only appear at higher order): they are purely nonlinear and describe a distortion of the spatial 
slices determined at this order by a constraint, quadratic in the scalar and vector variables. 

The main results of our analysis are as follows: (a) at leading order a purely Newtonian nonlinear 
energy current sources a frame-dragging gravitomagnetic vector potential, and (b) in the leading- 
order Newtonian regime and in the linear relativistic regime the two scalar metric potentials are the 
same, while the nonlinearity of general relativity makes them different. 

Possible applications of our formalism include the calculations of the vector potential [l],@] and the 
difference between the two scalar potentials from Newtonian N-body simulations, and the extension 
of Newtonian approximations used in structure formation studies, to include relativistic effects. 

PACS numbers: 98.80.-k; 98.80.Es; 95.35.+d; 95.36.-l-x 


I. INTRODUCTION 

The ACDM model [§, ;4j has emerged in the last few decades as the standard “concordance” model of cosmology || . 
Beyond photons and baryons, the main components of ACDM are cold dark matter (CDM), able to cluster and form 
structures, and the cosmological constant A, responsible for the observed acceleration of the Universe’s expansion. 

ACDM is based on Einstein general relativity (GR), and on the cosmological principle, i.e. a cosmological version of 
the Copernican principle [fi, 7j. This request for the Universe to be, on average , homogeneous and isotropic translates, 
in the language of spacetime, into assuming a Robertson-Walker metric. With this, Einstein equations reduce to the 
Friedmann equations; the solutions of these equations are the Friedmann-Lemaitre-Robertson-Walker (FLRW) models. 
Cosmic microwave background (CMB) anisotropy measurements [§i] have established that the Universe is, to a great 
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degree, spatially fiat, as confirmed by recent Planck [9] and baryon acoustic oscillations (BAO) data 0- In this 
standard scenario, small primordial inflationary perturbations on top of the FLRW background grow and produce the 
CMB fluctuations and the large-scale structure that we observe at low redshift. 

The theoretical tools that we use to study the growth of the large-scale structure are basically two: i) relativistic 
perturbation theory [ll|j EH is used to describe fluctuations in the early Universe, in the CMB, and in the matter 
density field on very large scales; ii) Newtonian methods, notably N-body simulations 0, are used to study the 
growth of structures in the nonlinear regime, at late times and small scales. Lagrangian perturbation theory (LPT) 
14, 15] (i.e. the Zel’dovich approximation or its second-order extension, 2LPT) is typically used to set up initial 
conditions for N-body simulations; LPT and other Newtonian approximations are also used to model nonlinear scales, 
e.g. BAO and CDM halos 

Observational cosmology has now reached an unprecedented precision, allowing stringent tests on models of the 
Universe. The tightest constraints come from the CMB 


22 ], 


supernovae 24 
scale structure 


A number of probes exist at low redshift, such as 
the local measurements of the expansion rate [24] and the three-dimensional mapping of the large- 
[25l|26|. In particular, rich clusters, lensing and redshift space distortion allow the measurement of the 


growth of clustering [27|,l28|. The standard ACDM model is very well supported by all these observations. However, 
a tension emerges in the framework of the base ACDM when parameter values measured from low-redsliift probes 
are compared with the values obtained from the CMB, as recently pointed out [29l-l32|. In particular, this tension 
shows up in a different growth of clustering in LSS at different epochs Q, as confirmed by the Planck release: “as in 
the 2013 analysis, the amplitude of the fluctuation spectrum_i_s found to be higher than inferred from some analyses 
of rich cluster counts and weak gravitational lensin g” (see [22] and Refs, therein). Measurements of redshift space 
distortion also suggest less clustering than the CMB [2E| 3JJ. Notably, adding a variable total neutrinos mass as extra 
parameter to the base ACDM model relieves the tension [30.[31.]. 

Given that the presence of a CDM component is widely accepted, while a cosmological constant has its own 
problems [H}, much work is gone into exploring alternatives to A. A number of possibilities have been proposed 
which, in essence, can be divided in three groups. One is to keep GR and modify the energy-momentum content, the 
dark sector in particular, replacing A with some form of dark energy (34] ]. A second alternative that has been widely 
explored in the last decade is to replace GR with a modified theory of gravity |35], 0. A third, more radical, option 
is that of abandoning the cosmological principle (37| , considering inho mog eneous models and the possibility that the 
observed acceleration of the Universe is the result of backreaction [H, [39j . either dynamical or optical [40]. 

Galaxy surveys are now aiming at the same 1% precision as CMB measurements: for instance, the detection of BAO 
by the Sloan Digital Sky Survey has recently allowed the first 1% level cosmological constraint by a galaxy survey [4l| . 


In addition, future surveys such as Euclid [42H44j and The Square Kilometer Array |45l - l47l ] will reach scales of the order 
of the Hubble horizon. It is, therefore, crucial that the theory used to make predictions and to interpret observations is 
developed with a matching accuracy. In particular, while exploring alternatives to the standard ACDM scenario is an 
interesting challenge, it seems timely to refine the theoretical modelling of ACDM, bridging the above mentioned gap 
between the Newtonian treatment of nonlinear small scales and the relativistic description of large scales. Ultimately, 
going beyond the Newtonian approximation in simulations of large-scale structure should be important in order to 
take into account causal, retardation and other GR effects that may be non-negligible for simulations that aim at 
~ 1% accuracy, in view of future surveys such as Euclid j42|, on scales of the order of the Hubble horizon. 

A first step in this direction would be to include GR corrections in the initial conditions for simulations, using a 
dictionary based on first-order perturbation theory (48], cf. j4f|; also, first-order GR effects on horizon scales must be 
taken into account in interpreting bias and non-Gaussianity |50[. However, while the Poisson equation in Newtonian 
gravity establishes a linear relation between the gravitational potential and the matter density field, the intrinsic 
nonlinearity of GR unavoidably generates new effects, even when perturbations are small. For instance, an initially 
Gaussian curvature inflationary perturbation translates into an effective non-Gaussianity of the density field in the 
matter era (5ll - f54l ]. a nonlinear effect that should be taken into account in initial conditions for simulations. But the 
Poisson equation embodies action-at-tlie-distance; i.e., it is the mathematical representation of the acausal nature of 
Newtonian gravity. Thus the ultimate step forward would be that of investigating the effects of relativistic nonlinearity 
in cosmological structure formation, including GR corrections in the evolution of the matter density field in N-body 
simulations, as well as in approximate treatments of the nonlinear regime 0SIIM. 

Aim of this paper is to present a new nonlinear relativistic post-Friedmann (PF) formalism which, in essence, 
is a generalization to cosmology of the post-Minkowski (weak-field) approximation, married with the fundamental 
assumption of the post-Newtonian (PN) approximation that velocities are small. Our goal is a relativistic framework 
valid on all scales, and including the full nonlinearity of Newtonian gravity at small scales. In this framework - 
assuming a flat ACDM background and a fluid description of matter - the exact equations of Newtonian cosmology 
[fill [62 ] appear as a consistent approximation of the full set of Einstein equations, determining leading-order terms 
(which we call OPF) in the metric. GR corrections (which we call 1PF) appear next, and are of two types: terms 
quadratic in the Newtonian variables (we may call these proper PN terms), and proper linear GR terms. Once 
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linearised, the nonlinear approximate equations we obtain for a set of appropriately resummed variables are those for 
the scalar and vector sectors of first-order relativistic perturbation theory [63|. Tensor modes are purely nonlinear 
and nondynamical; i.e., they do not describe propagating gravitational waves, but rather a distortion of the spatial 
slices. 

Let us clarify how our PF formalism differs from the traditional PN approach 64] and its various applications in 
cosmology. In a contemporary perspective [HJ, the correct derivation of the post-Newtonian approximation on the 
flat background spacetime follows consistently from the post-Minkowski approximation once the v/c <C 1 assumption 
is made, which also implies to neglect time derivatives with respect to space derivatives 1 . However, in cosmology the 
background is a FLRW solution, in our case the flat ACDM model, and what we can assume to be small are peculiar 
velocities, not the change with time of the physical distance between two arbitrary observers. To illustrate the point, 
adopting a Newtonian perspective, in the absolute space of Newtonian cosmology one uses the background comoving 
coordinates x as the Eulerian grid, and the physical pos ition of a fluid element is r = ax. Then r = Hr + v, where 
v = ox is the physical (or proper) peculiar velocity |ll|, representing the deviation from the Hubble flow; x is the 
peculiar velocity with respect to the comoving grid [62(. Then, should we assume that |i’| -C c, we would end up with 
an approximation only valid at small scales, well below the Hubble horizon, |r| <C off -1 . In addition, traditionally the 
PN formalism has been developed to study GR corrections to the orbits of isolated objects. Thus, the focus is on the 
equations of motion, rather than on a consistent approximation of the full set of Einstein equations. In developing his 
PN treatment of relativistic hydrodynamics, having in mind applications to relativistic stars, Chandrasekhar [67| also 
focused on the equations of motion for the fluid, as derived from the conservation equations. In both cases, the metric 
and other variables are expanded 2 in powers of 1 /c, then the expansion is applied to the equations of motion. As it is 
well known, a correction to the time-time component of the Minkowski metric is all is needed to obtain the Newtonian 
equations of motion. All other corrections are then considered post-Newtonian. We will discuss this point in detail 
in Sec. IVIB1 where, following [69], we call this approach focused on the equations of motion ’’passive”. In cosmology, 
however, we find desirable to consider an ’’active” approach, where the space-space and time-time components of the 
metric are equally weighted, leading directly to a consistent treatment, order by order in the expansion parameter 
1/c, of the full set of Einstein equations 3 . In view of applications, a useful byproduct of this active approach is that 
the equal weighing of the different metric components is what is needed to obtain the correct photon trajectories, i.e. 
the null geodesics that are at the base of the causal structure of the spacetime. 

Various early and more recent works have applied the traditional PN exp ansion in cosmology: in HHzIl PN 
equation of motion for particles in the expanding Universe are derived; in 74] the authors clarify the role of the electric 
and magnetic part of the Weyl tensor in the Newtonian approximation; in (75|_and [76| the PN analysis is given in 
Lagrangian coordinates using, respectively, the 3+1 and 1+3 framework; in [771 ] the authors derive a complete set 
of field and hydrodynamic equations in Newtonian-like forms. Other remarkable works following a PN method are 
those of Szekeres [78, 79). Finally, (8(| follows a ’’hybrid approximation scheme”, somehow closer in spirit to our PF 
approach, a mix between standard cosmological perturbation theory and PN approximation, also obtaining equations 
for the generation of gravitational waves. 

All these works apply the standard iterative approach of the PN expansion. Here instead we focus on defining a 
set of resummed 1PF variables, and on deriving the set of nonlinear evolution and constraint equations they satisfy. 
Specifically, assuming the Poisson (or conformal-Newtonian [HI, HI, |63j, HU ) gauge, we expand the metric in powers 
of 1/c and in the equations we retain all scalar terms up to order 1/c 4 (and all vector terms up to 1/c 5 ), instead 
of peeling-off the different orders, cf. 79]. This leads to a final set of equations that differs from those of the 

works mentioned above. Partly following the PN tradition, the OPF and 1PF orders, respectively refer to terms 1/c 2 
and 1/c 4 , relevant for scalar potentials. The OPF terms are Newtonian, the 1PF terms contain the GR corrections. 
Finally, we stress that, in order to obtain a system of equations with a well-posed Cauchy problem, we should consider 
terms of the next 2PF order, i.e. 1/c 6 [78j, 79|. 

In summary, the main goal of this paper is to present a set of nonlinear resummed equations up to 1PF order which 
retain in full the nonlinearity of Newtonian theory on small scales and all linear relativistic perturbation theory on 
large scales. The 1PF scheme is, therefore, capable of describing, in a unified framework and at the relevant leading 
orders, the evolution of large-scale structure on all scales of cosmological interest. Our two main results simply follow 


1 The contemporary derivation of the post-Newtonian approximati on |65| should be contrasted with the traditional implementation (see e.g. 
fan which suffers from various ambiguities and inconsistencies |65l l66H . By and large, cosmological applications of the post-Newtonian 
approximation follow the traditional approach (see below). 

2 Formally, the expansion is with respect to a small dimensionless parameter representing the ratio v /c; by dimensional analysis the ratio 
of the gravitational field and e 2 is also dimensionless and can be expressed in terms of the same parameter. In this way the parameter, 
which in practice is the inverse of the speed of light 1/c, indicates the relativistic weight of each term in the expansion [68 ]. 

3 The approximate conservation equations then consistently follow from the contracted Bianchi identities, as in the exact theory. 
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from the analysis of our nonlinear equations: at leading order a purely Newtonian nonlinear energy current sources a 
frame-dragging gravitomagnetic vector potential; in the leading-order Newtonian regime and in the linear relativistic 
regime the two scalar metric potentials are the same, while the nonlinear 1PF equations for the resummed scalar 
potentials imply that the nonlinearity of GR makes them different. 

The paper is organized as follows. After defining the various metric terms in Sec. UH we obtain the stress energy 
tensor in Sec. Mil the field equations in Sec. IIVI and mass and momentum conservation equations in Section [Vl We 
then consider our equations in two opposite limits: the Newtonian regime on small scales, neglecting 0(l/c 4 ) terms 
(see Sec. I VII) and the linear regime on large scales through the linearization of the equations (see Sec. IVIII) . In Sec. 
IVIII1 we first define suitable resummed variables, then we obtain a consistent set of nonlinear equations describing 
their evolution. In Sec. m we draw our main conclusions. Finally, in Appendix [AJ we apply the PF expansion to 
the Riemann and Ricci tensors. 


II. NEWTONIAN AND POST-FRIEDMANN METRIC VARIABLES 

We consider a homogeneous and isotropic FLRW flat background where two kinds of perturbation terms are added to 
the metric, representing two different levels of accuracy. The first will give the Newtonian regime, i.e. an approximate 
solution of Einstein equations such that the dynamics is described by the exact nonlinear equations of Newtonian 
cosmology for pressureless matter. The second will give the relativistic corrections, adding terms to the Newtonian 
equations. In our post-Friedmann framework the goal is to calculate these relativistic terms; in this spirit, we shall 
refer to the Newtonian approximation as the OPF (or leading) order, and to the first relativistic corrections as 1PF 
order. 

The expanding parameter for the inhomogeneous perturbations is formally given by 1/c, so that the components 
of the metric tensor in the line element 


can be written as 


— c 2 dr 2 = ds 2 = g tlly dx l ‘dx v 


g oo 

got 

gij 


1 — — 5 — + -j{2Un — 4t/p) 


+ 0 


= a 


a tdN a tdP 

2V n 


O 


1 + — 3 - + -^(2V n + 4Vp) ) 5ij + -jh ZJ 


+ 0 


( 2 . 1 ) 


(2.2a) 

(2.2b) 


(2.2c) 


where r is proper time and a(t) is the scale factor of the FLRW background. Note that we assume that both the scale 
factor and the metric are dimensionless, while the coordinates have dimension of a length. Greek indices take the 
values 0,1, 2,3 and refer to spacetime coordinates, Latin indices refer to the spatial coordinates. In particular in our 
post-Friedmann scheme, for a proper “powers of c” order counting, it is important to note that the time coordinate is 
x° = ct. Here the Kronecker Sij represents the metric on the flat spatial slices of the background; the spatial Cartesian 
coordinates are understood as an Eulerian system of reference, cf. [821,183}. 

Our metric is a generalization to cosmology of Chandrasekhar’s metric for post-Newtonian hydrodynamics [67} (cf. 
s for a cosmological application). However, it is important to remark the difference of our post-Friedmann approach 
from the standard post-Newtonian one |6J, 1651, |8J, |85[ . In the latter, the focus is on the equation of motion for matter, 
hence only the leading order g 00 metric perturbation is Newtonian, while the is post-Newtonian. In our approach 
the focus will instead be on the complete set of Einstein equations. We aim at a self-consistent approximate set of 
equations at each order, similarly to the post-Minkowski (weak field) approximation [65l. [86]. Then, consistency of 
the Einstein equations dictate that the gij and goo metric perturbations must be of the same order (see Sec. IVII for 
the Newtonian approximation). Anticipating these results, the indices N and P label Newtonian and post-Friedmann 
quantities. The N quantities are the only relevant one at the OPF order of approximation of Einstein equations, i.e. 
in the Newtonian regime. The P quantities appear at the 1PF order. 

We assume the Poisson (or conformal Newtonian) gauge [ill, d, HI, IHI: H3] , so that the three-vectors Bf and Bf 
are divergenceless, Bf' 1 = 0 and Bf’ 1 = 0, and hij is transverse and tracefree (TT); i.e., it represents pure tensor 
modes h\ = h r ,- 1 = 0. Commas in front of indices have the standard meaning of partial derivatives. 

Having completely fixed the gauge leaves us with six degrees of freedom at each order: the two scalars Un and 
Vjv at leading order, and Up and Vp at 1PF order; the two independent components of R) v (Bf at 1PF); the two 
independent components of hij. However, the latter only appears in the equations at 1PF order. At leading order, 
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Einstein equations impose C/.y = V)v, i.e. a single scalar gravitational potential in the Newtonian regime, as expected. 
At leading order, is determined by the (vector part of the) Newtonian energy current, and cannot be set to zero. 
It is not dynamical; i.e, it doesn’t contribute to matter motion. However it appears in the metric and does affect null 
geodesics and observables. It can, therefore, be extracted from Newtonian N-body simulations 00 and it contributes 
to lensing 


Note that the 1PF corrections consist of quadratic combinations of Newtonian quantities and intrinsic 1PF variables. 
We have also included tensorial TT modes but, at 1PF order, they cannot be interpreted as gravitational waves; they 
satisfy a constraint equation rather than an evolution equation (see Sec. IIV C 31 c.f. [83jl. 


III. MATTER VARIABLES 


Having defined the metric variables and their weight with respect to the expansion parameter c 1 , we now look at 
the matter quantities. The dimensionless 4-velocity is defined as 


and it satisfies the usual relation 


dx» 

cdr 


g^u^u v = - 1 , 


(3.1) 


(3.2) 


i.e. is a unitary timelike vector field. In Newtonian cosmology one uses the background comoving coordinates x as 
the Eulerian grid, the physical position of a fluid element is r = ax, and v = ax and x respectively are is the physical 
(or proper) peculiar velocity |Tl| (the deviation from the Hubble flow) and the peculiar velocity with respect to the 
comoving grid [62i ]. With this in mind, it is then natural to define the physical peculiar velocity as v' := adx l /dt. 


Therefore, we also define Vi := 5ijV J . 
Then, 


dx l 


dx l dt v l 


cdr cdt dr ca ’ 

and, using dSUD and keeping terms up to order c 4 , the 4-velocity components are 

, 1 v i 


u " =- u°, 

c a 

1 


1 


1 


u — 1 H—( Um + — v I 4—j 


u 0 = -1 + 


--B™ + ViU N + 2viV N + -ViV 


)-U 2 N + 2 Up + v 2 V n + lv 2 U N + V - 


1 


U N - lv 2 ) + 


2 Up - - \v 2 U n - v 2 V n - 

Z Z o 


(3.3) 

(3.4a) 

(3.4b) 

(3.4c) 

(3.4d) 


We consider a Universe filled by cold dark matter (CDM), described by a single pressureless (dust) component with 
energy-momentum tensor 


= c 2 pu^Uv, 

where p is the mass density. 

In the Poisson gauge the components and trace of T 1 '^ then are 

1 


T° 0 = - c 2 p-pv 2 -^p[2(U N + V N )v 2 -B?v l + v A ] 

T°i = cpavi + -pa{vi[v 2 + 2(U N + V N )]-B^} , 
c 

T l 0 = — c—pv l — --pv 2 v l , 
a c a 

T*j = pv'vj + \p {uSV + 2 (E/jv + V N )} - v l Bf \ , 


T\ = T= -pc 


(3.5) 

(3.6a) 

(3.6b) 

(3.6c) 

(3.6d) 

(3.6e) 
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All these quantities are written in order to explicitly show the different contributions in powers of c, where in each 
expression the first term beyond the background term (if present) represents the leading Newtonian order, the second 
the first 1PF correction, etc. Note that there is no approximation in the trace: indeed, the mass density p plays the 
role of a fundamental exact quantity that is not expanded into contributions at different orders. In the following it 
will be useful to use the density contrast 6 , defined as usual: 6 := (p — p)/p, where p denotes the background matter 
density. 


IV. EINSTEIN EQUATIONS 
A. Expansion of Einstein equations at 1PF order 

We now consider Einstein field equations for the metric (12.211 . including the cosmological constant A: 


= A"„ - lR6» = - AS fl „ 


(4.1) 


Expanding in powers of 1/c, in all equations we retain the first two terms of the expansion. We then obtain the 
following equations, where the dot denotes partial differentiation with respect to coordinate time t. 


Time-time component 
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Time-space component 
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Finally, neglecting in Eqs. m and (14.81) all terms representing inhomogeneities, we obtain the background equations: 



1 8t tG 
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These are recast into the standard Friedmann and Raychaudhuri equations for the flat FLRW background 

1 rr2 1 8 t rG_ A 

? H + 

H + H 2 


l r 

72 


1 47rG 


3 
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c 2 3 


"P + 7T j 


(4.6a) 

(4.6b) 


after substituting the Hubble expansion scalar H = a/a. 

B. 1PF equations for the inhomogeneities 

We now subtract the background parts (14.511 from Eqs. (14.21) and (14.31) . in order to obtain equations for the inho¬ 
mogeneous quantities. The time-time component of the field equations then gives a generalized Poisson equation: 




1 47rG . 1 47rG .. 2 
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(4.7) 


Note that the cosmological constant disappears from the equations above; it only directly contributes to the back¬ 
ground dynamics, i.e. Eqs. (14.61) . Thus perturbations are only affected by A through their coupling with the Hubble 
(background) expansion. 

The trace of the space-space component (TOll gives 
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while the trace-free part is 
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C. Scalar, vector and tensors parts 


1. Scalar equations 

It is useful to recast the previous equations in order to isolate, in the linear part of the equations, the scalar, vector 
and tensor contributions. For the scalar sector let us apply the divergence operator on Eq. (TOt : 




2 V 2 | ^Up + Vpj + (v N Ui) ’* + 2-^ ( U N Ut)’ 1 + 2 (Vn^Vn) ' 
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Multiplying Eq. (14.101) by H/c, and applying V 2 on Eq. (14.71) . we obtain the following constraint equation: 


(4.10) 


1 9 9 1 

—V 2 V 2 Vv - -j 


V 2 V 2 (V£) - -V 2 (V N ,iV N ■<) - 2 V 2 V 2 Fp 


= — 47 rG(a 3 p)] — -V 2 S 


iv 2 ((l + <5)^ 2 )-3-(u i (l + 5))’ i 

a a 


(4.11) 


Now, we can obtain a second scalar constraint by applying the operator djd 1 on both sides of Eq. (|4.9p . This then 
gives: 

12 V 2 V 2 (V n -U n ) + ^{ Vv 2 (Vp - Up) + (u Nti u N ’ j ) ’’ - (v N>i v N < j ) ' l + 2 (u N>i v N ’ 2 ) 


c 2 3 

1 . 


Jv 2 (u Ntk U N ' k ) + iv 2 (v N , h V N ’ fe ) - ^V 2 (u Ntk V N ' k ) - 2 |Vjv(Viv - U N )/ 


J + 3 V 2 [V n V 2 {V n -U n )] 


8 ttG 


—a 2 p 


(1 + 5) ( viv 3 - — v 2 


J ,3 


(4.12) 


These equations will be useful in Sec. IVIIII 


2. Vector equations 

The vectorial part of equation (14.41) can be found using the curl operator: 

-4v x v 2 £f + 4 V x [-^Bf + b? k (v N - u N y k - b£ ,(u n + v N y k + 2 v 2 b»{v n - u N ) + b?v 2 v n 

C° C° ’ ’ 

167rGa 2 „ . . 167rGa 2 „ , , r 9 nK1 , 

= - 3 —V x ( pvt ) H-g— V x (p {uj [u + 2(Un + Viv)] — B^ }) . 


+2 B?V n , 


(4.13) 


Alternatively, we can also obtain a constraint equation for the vector part by applying the operators V 2 to Eq. m 
and dj to Eq. (14.101) . finding 


1 V72 V72 toN i ^ 


9A r v72l 


{v 2 v 2 £f-v 2 

>k (V N - 

+ V 2 


(Viv + GAr)’ fc 

- 2V 2 

+ 2 

B? >k (V N - Giv)’ J ' fe ' 

+ 

V 2 Sf , fc (Vat - 


72 D N , 




+3 (BfV 2 V N + 167 rGa 2 p (<WBf) . - 167 rGa 2 p V 2 [(1 + 6) Bf] 167 rGa 2 p {[^(1 + 5)]j - V 2 [v t { 1 + 5)]} 

,k 


+?i 4 “ v 

+167rGa 2 p 


( 'VnUn - (VwUjV'kJ + 2-V 2 (UnUn,i) - 2- (t/jvf/jv.fc)’* + 2V 2 (VnjVn^ — 2 (Vn^Vn'J 


ijil + S) (v 2 + 2 U n + 2V n )’ 


- 16nGa 2 p V 2 


z (1 + <5) (^ 2 + 2{7 jv + 2Viv) 


(4.14) 
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These vectorial equations also depend on scalar quantities, because of nonlinearity; however, the divergence of these 
equations would give zero, as it should. 

An evolution equation for B, can be obtained by applying the operator djV 2 on Eq. (14.911 and subtracting the 
equation obtained applying d t on Eq. (14.121) . This procedure leads to 


1 • 


^V 2 V 2 [ ^Bf + -Bf ) =__<iV 


JV 


2 1 a 1 


1 1 


(tWV) . -(u N , k u N - 


- v- 


(v N , 


V 


(v N , k v N A 


+v- 


(u N ,iV N + v 2 - 2 (B A r, fc y v J ) ’* 


— 2V 2 


8 ttG 

V ' 


(1 + 5)(wit;* - -dfc 2 ) 


- ,fc 


(Vat — Un)’^Vn 
,k 


+ 2 


J j 


(V N - Un)’{Vn 




(1 + S)(v k v J - -5 3 k v 2 ) 


(4.15) 


This shows that an evolution term for BA only appears at order 1/c 4 . 


3. Tensor equations 


In order to isolate the TT part of the metric hij we define the following nonlinear quantities: 


-V = U Nti U N ’ 3 - V N ,iV N ’ 3 + U Nti V N ’ 3 + U N ’ 3 V N ,i - 2V n {V n - U N ) f 


s 3 


- ~^UN,kU N ' k + -VN,kV N ' k — -UN,kV N ’ k + -VatV 2 (Vat — Un) 


<V = (1 + 5) 


ViV j - 


With these definitions, Ea. (14.15l) becomes 


-V 2 V 2 


2 a 


-B 


N 


1 a 


Bf I = 


(v-U)., • s^/p( 




Finally, using Eqs. (14.91) . (14.121) and (14.181) . we obtain the following constraint equation for hij-. 


i-V 2 V 2 V 2 /y = \ 


-A l £% - V 2 A l k k S 3 + 2 V 2 A k i + 2X7 2 A k M 6 l3 - 2V 2 V 2 A/ 


+87rGa 2 p - ^ 2 S k k S 3 + 2V 2 S? M 6 lj + 2V 2 S k ’ 3 - 2 V 2 V 2 S) 


(4.16) 

(4.17) 


(4.18) 


(4.19) 


V. CONSERVATION EQUATIONS 

The field equations in GR are constructed in order to imply the conservation equations through the contracted 
Bianchi identities (89): 


T 


V 


= o. 


(5.1) 


From this, considering the time component and keeping all terms up to 1/c 2 order, we obtain the energy 
conservation equation 



OV),i 


(■ a 4 pv 2 )' 


3 VnP + 


( pv 2 v i ), i 


pv 


(3 Vjv — C7jv),» 


= 0 


(5.2) 


Setting v l = Un = Vn = 0, and p = p, this equation reduces to the background continuity equation for cold dark 
matter, i.e. p = —3Hp. Subtracting this equation from (15.21) and defining the total, or convective, derivative for any 
quantity Q as 


dQ _ a v l Q,i 

dt a 


(5.3) 
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Eq. m becomes 


dS 

dt 


+^(5+i)+i 

a c A 



CL o 
-V 2 + 

a 


d(v 2 ) 3 dV N 

dt dt 


- v -u 


NA 


= 0 . 


(5.4) 


Note that we have used the Newtonian part of this equation to simplify the 1PF order part. The space part of (15.11) 
gives the momentum conservation equation: 


(a 4 /w;)' pU N}i ( v j pvi) tj 


A2,V N -U N ) tj pUp i {pvivH 2 )^ ( a 4 pv 2 Vi)' 

pviv j ---2- ’- + pvi(3V N -U N ) H-— H- j - 


[2Apvi(V N + U N )]' [pv i v i (V N + U N )] i . 2 {V N + U N ) ti (a 4 Bfp)' (BFpv*) 


+ 2 


B 


N 


pv 


H-— pv J > = 0 . 


(5.5) 


Note that Eqs. m and (15.51) are equivalent to Eqs. (57) and (58) in & as far as the conservation equations are 
concerned, there is no difference between the standard post-Newtonian approach used in [77} and our post-Friedmann 
approach. The difference becomes relevant for the consistency of the full set of Einstein equations; we discuss this 
point in Sec. IV1 Bl Simplifying this equation by using the Newtonian part of Eq. (15.41) and the background continuity 
equation for p , we derive the 1PF Euler equation: 


dvi a Un,i ] 
dt + a 1 a + c 

1 d Ns Bf iV i 

a dt a 


Vi —-(Un + 2 Vn) + -Un,z(Un + V)v)- v 2 Vn,i - Upi + -v^Unj - -v 2 Vi 

dt a a a a a 


= 0 . 


Finally, using the Newtonian part of Eq. m, the continuity equation can be recast in the following way: 

= 0 . 


dS vl i,c: ^. 1 


(£+!)( ~v j U N j - -v 2 + 

' a a dt 


(5.6) 


(5.7) 


The last two equations are the hydrodynamic equations of motion of the dust component in the post-Friedmann 
approximation, at 1PF order. Let us stress that, in Eqs. (E3 and (15.61) . the post-Friedmann corrections contain 
nonlinear coupling terms between metric variables and matter variables, p and u,;. This is the reason why the vector 
potential Bi cannot be decoupled from the scalar modes in the equations of motion. 

In the next two sections we will check the consistency of our approach with our goal of obtaining a set of equation 
valid at all scales. To this end we will first consider the leading order, the Newtonian regime, which is a good 
approximation on scales much smaller than the Hubble radius, then we will linearize our equation to check consistency 
with linear relativistic perturbation theory in the Poisson gauge fTTl. Ktl . which is a good approximation on large scales. 


VI. LEADING ORDER: THE NEWTONIAN REGIME 

A. Newtonian dynamics from consistency of Einstein equations, with a bonus 

Retaining the leading order-terms in the c -1 expansion, i.e. the OPF order, we recover the equations of Newtonian 
cosmology and we obtain the corresponding spacetime metric. We call this the Newtonian regime: as we are going to 
see, the dynamics is purely Newtonian, yet we have a spacetime metric that is a well-defined approximate solution of 
Einstein equations. 

We obtain the Newtonian continuity and Euler equation from the hydrodynamic equations (15.41) and (15.61) : 


<5 + 

Vi + 


v l Si 


H- -(6 + 1 ) = 0 , 


v°Vi,j a 1 
- 1 — Vi — —Un , 


( 6 . 1 ) 

( 6 . 2 ) 
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At leading order, from Einstein equations (14.71) . (14.41) . (14.81) . (14.91) . we obtain 


G u 


■ A = 


8 t tG 


T °o -> 44 v 2 ^v = - 


4ttG 


pS 


G°. = 


trace of G i j + A S i 3 

8 t tG 


~ A + 2^^ + -VV.i 

2 ar a z a 


\^V 2 (V n ~U n ) = 0 


8 ttG 


p(l + d)uj , 


1 1 


traceless part of G % j + AS 1 3 = ——TV ->• — — [(Vjv - Gjv)^’- 3 - -V 2 (Vat - C/jv)^] = 0 . 


c er 


(6.3a) 

(6.3b) 

(6.3c) 

(6.3d) 


Equation (|6.3al) is the Poisson equation, from which the Newtonian character of the spatial metric potential Vjv is 
apparent as long as we identify it with the Newtonian gravitational potential generated by the matter field, Vjv = —< f>N■ 
As noted by Bertschinger [86|, it is an interesting fact that at leading order only g-ij contributes to G° 0 . On the other 
hand, as is well known, it is the time-time metric potential Um that tells matter how to move in the Euler equation (16.21) 
(or to particles in the geodesic equation; see e.g. [69}). Equations (I6.3cl) and (I6.3dl) show that Eqs. (14.81) (14.91) reduce, 
at leading order, to constraint equations: they give the consistency relations between the scalar metric potentials f/y 
and Vjv in a GR context at this order. Modulo residual gauge modes [86j] that leave (16.3d) and (I6.3dl) invariant, we 
must have Un = Vv = —4>n- In summary, Einstein equations reduce, at leading order, to the standard equations 
of Newtonian cosmology. The metric tensor generated from a self-consistent expansion of the full set of Einstein 
equations at leading order is the cosmological version of the weak field metric, with the FLRW metric replacing 
Minkowski as background. 

Our metric (12.21) . however, also contains the frame dragging vector potential B t and the TT part h l3 . This TT 
part only appears at 1PF order in G* and so it is irrelevant at leading order, consistently with the fact that, in this 
Newtonian regime, we only retain the Vy term in gij. On the other hand, the leading order of the G 3 equation is 
Eq. (I6.3b[) : thus, it can not be neglected. This is a new equation that arises in the relativistic context and determines 
the leading-order frame dragging potential from the purely Newtonian term on the right hand side, the energy 
current that is determined by the other equations. Therefore, even in the Newtonian regime, the frame dragging term 
cannot be set to zero. Some authors consider this term of higher order; however Eq. (I6.3bl) tell us that this is 
inconsistent with the other Newtonian equations. Indeed, taking the divergence of (I6.3bl) (i.e. its scalar part) and 
using the Poisson equation implies the continuity equation, while taking the curl (i.e. the vector part) shows that the 
curl of Rj v , a gravitomagnetic field [86}, is sourced by the curl of the energy current pVi, and there is no reason why the 
transverse part of this current should vanish. Alternatively [86], one can split vectors into longitudinal and transverse 
parts. Then one can see from (I6.3bl) that a purely Newtonian transverse energy current sources BjA. Thus, setting 
B^ = 0 does imply, even in the Newtonian regime of dynamics, a purely longitudinal energy current, i.e. imposing 
the condition = 0 does imply an artificial constraint on the general Newtonian dynamics. While in the linear 
regime vectors modes can be separated from scalar modes, in the nonlinear regime of Newtonian dynamics there is no 
reason why the energy current should be longitudinal, and indeed it is well known that transverse vector modes are 
generated by nonlinearity, and in particular vorticity is generated after shell crossing [90]. The nonlinearly generated 
gravitomagnetic potential B™ can actually be extracted from standard Newtonian N-body simulations [2, @], and 
its contribution to lensing has been computed in [88]. The bottom line is that in the Newtonian regime of GR the 
gravitational field is not purely scalar, it does contain a gravitomagnetic vector part urn- 

As we are going to show in the next section, the above considerations on the frame dragging Bi are entirely 
consistent with linear perturbation theory. In the perturbative regime one naturally splits vectors into longitudinal 
and transverse parts and finds that Bi decays [78H8CI |92| . [93]]. Nonetheless, in the nonlinear Newtonian regime the 
transverse part of the energy current sources Rj v |jJ, ]2|] . At this order, however, R) v is a nondynamical variable, rather 
it is determined nonlocally through the acausal equation (16.3bl) . Furthermore, a calculation of the magnetic Weyl 
tensor shows that, at leading order in powers of c~ , 


Hij 2c 3 


B 


'fj,,v(i £ j) ^ + 2u m ((7at + Viv ),v{i £ j) 


fiiy 


(6.4) 


so that its linear part, i.e. the first term in the square bracket, is precisely the curl of B™. As discussed in detail 
by Kofman and Pogosyan [74[, although the magnetic Weyl tensor is zero in Newtonian theory (i.e. starting from a 
scalar theory of gravity), it cannot be neglected in deriving the “covariant equations” IH in the Newtonian regime 
from relativistic theory. Nonetheless Hij is not an independent dynamical variable in the Newtonian regime, rather 
Eq. (16.41) shows that it is locally determined by other variables. 
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B. Passive and active approach 

Starting with Einstein |89j, the usual way to consider the Newtonian limit in GR consists in demanding that the 
timelike geodesic equation agrees with the Newtonian equation of motion for a particle j64[ f69| . This naturally sets 
<7oo = —(1 — 2Un/c 2 ), identifying — Un with the gravitational potential cj) N . For a fluid, as we can see from (16.11) 
and (16.21) . the continuity and the Euler equations do not contain Vn, so that the latter could be interpreted as a 
post-Newtonian variable. This way of thinking considers the “passive” aspect of gravitation, namely the response of 
matter to gravity [H]: one wants to determine the equation of motion of a particle in a given gravitational potential. 

In addition, the Newtonian limit is again considered using the fiel d eq uations with matter, in order to determine the 
gravitational coupling constant by recovering the Poisson equation [89|. Starting from the following form of Einstein 
equations 

= k(T^ v - ^TSftv) + KSfjtv (6.5) 

one finds that the time-time component only contains goo , so that, at 0(c~ 2 ), we obtain the Poisson equation for Un 


1 9 r r 47tG 

V U N =- „ r, PS 


3 a 2 


3 c 2 


( 6 . 6 ) 


by identifying n = 8irG/c 4 . Again, one would be tempted to conclude that only Un is needed at leading Newtonian 
order. 

However, in our approach, we are interested in the “active” aspect of gravitation, i.e. the generation of gravity by 
the matter distribution. To this end, we need to determine the metric using the field equation in a fully consistent 
manner, even in the Newtonian regime. In particular, in cosmology we deal with a self-gravitating fluid, so that the 
complete set of field equations is as important as the equations of motion. 

In the passive Newtonian approach of the standard post-Newtonian formalism the geodetic equation determines 
the order of the metric variables, Un as Newtonian and Vn as the post-Newtonian correction. In our post-Friedmann 
scheme it is the leading order of the field equations, derived from the complete set of Einstein equation, that establishes 
the order of the metric perturbations, giving that both Un and Vn have a Newtonian character. This emerges from 
(I6.3al) . (16.3cl) and (I6.3dl) and naturally also from the spatial components of (16.51) . whose trace is 

— — V 2 (Un~4:Vn) = —j— P& ■ (6-7) 

&■ Or &■ 


This clearly shows that neglecting Vn with respect to Un would not yield the Poisson equation. 

In summary, if we want to obtain the equations of Newtonian cosmology from GR, in the active approach and 
consistently considering all components of Einstein equations, we do need both scalar potentials Un and Vn as OPF 
terms in the metric, with the result 


Un = Vn = — 4>n ■ 


( 6 . 8 ) 


VII. LINEARIZATION: RECOVERING FIRST-ORDER PERTURBATION THEORY 

Standard relativistic perturbation theory [lj, UH, H3, HU HI, HU , where density, velocity and metric variables are all 
assumed to be small, is a fundamental tool in theoretical cosmology. It provides the framework to develop predictions 
from inflation for the matter and metric fluctuations in the early Universe and to work out their imprint on the CMB 
and the matter fluctuations at the beginning of the matter dominated era [53] • 

Our post-Friedmann scheme has been developed to generalize the equations of nonlinear Newtonian cosmology to 
include relativistic corrections. It would, however, be a great bonus if, in linearizing our equations, we could recover 
standard first-order relativistic perturbation theory. Given that the observed Universe is remarkably well described 
by perturbation theory at large scales, up to the Hubble horizon and larger, recovering at least the first order would 
imply that our scheme is not only valid at small nonlinear Newtonian scales, as shown in the previous sections, but 
also at the largest scales of interests in cosmology. Note that it is not a priory obvious that this recovery is possible; 
the aim of this section is to explicitly show that, by linearizing our equations and by defining resummed metric 
variables that contain a Newtonian part and a relativistic 1PF correction, we obtain first-order perturbation theory. 
This resummation will also be at the base of the following sections. 
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In first-order perturbation theory, it is standard to decouple the scalar, vector and tensor modes [II], OH (63], |H, 
!H, It is, therefore, convenient to separately look at the scalar and vector type equations of Sec. II V Cl las already 
noted, at 1PF order the tensorial modes are not dynamical, i.e. we can’t recover first-order gravitational waves). 
In particular, the first-order velocity perturbation v l can also be split into a scalar and vector (solenoidal) part, 
v l _ v || ,* _|_ v i ± ) where v]_ ^ = 0. 

Starting from the scalar sector and considering only the linear terms, we obtain from equations (14.71) . (14.811 . (14.101) . 
and (14.121) 


1 


a • 

( a\ 

-v N + - 

a 

\ a J 


} V 2 Vp 

J 3a 2 


1 4t tG 


a 2 p5 


[V 2 (U n - V N )] + — < 2V 2 (Up - V P ) + 3a 2 


(Un + 3Vn) + 2 — Un + ( — ) Un + Vn 
a V /a V a 


-iv 2 ( a -u 


N + Gy ) + 4 (2-V 2 [/p + 2V 2 Vp ) = ^AirGapd , 


1 r V 2 V 2 (V) v - J7jv) + 12V 2 V 2 (Vp - Up) = 0 , 


(7.1a) 

(7.1b) 

(7.1c) 

(7-ld) 


where, in Eq. (17.1cl) . we have defined the spatial velocity divergence as 0 = v l i = V 2 r^ . The continuity equation and 
the divergence of the Euler equation from (E3 and m , after linearization, are 


S H-1—j Viv — 0 , 

a c z 

■ , d a V 2 U n 2 V 2 ?7p 

t/H— v -o- 

a a c z a 


= 0 , 


(7.2) 

(7.3) 


where the convective term in the total time derivatives is neglected and partial and total derivatives coincide. Finally, 
defining the following resummed scalar metric variables 


< j>p := — {Un + -zUp), 


"0 p '■= + ~2 Vp) 


the previous equations become 


V 2 i/j p - -a 2 


-4>p+ - </>p 


= A-KGpa 2 8 , 


—V 2 (i/’p — 4>p) + —a 2 ~(4>p + 3ipp) + 2— (frp + ( — ) 4>p + ipP 


= 0 , 


V” 0p + ippJ = —4irGapO , 

4 v 2 V 2 ( 0 p-V'p) = 0 , 
c 2 

i e 3 • 

<5-1-2 t/’p = 0 , 

a c 

0 + -9 + iv 2 </>p = 0 . 
a a 


(7.4) 

(7.5) 

(7.6a) 

(7.6b) 

(7.6c) 

(7.6d) 

(7.6e) 

(7.6f) 


Note that, from Eq. (17.6dl) . ip-p = cf>p\ this follows from the perfect fluid matter model we are considering (i.e. null 
anisotropic stress, see e.g. [93)• The equations above coincide with those for scalar fluctuations of standard 

perturbation theory in the Newtonian gauge; see e.g. [Ill, [6^, HU . The conclusion of this analysis is important: the 
1PF approximation, in the linear regime, completely includes first-order scalar relativistic perturbation theory, so that 
all the linear GR terms are just 1PF. On the other hand, the 1PF order does not involve just linear terms, indeed 
there are contributions of second and higher order in the standard perturbative expansion. 
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This result is due to our active approach, where Un and Vn are both Newtonian variables, to the derivation of 1PF 
equations that retain the first two orders in the 1/c expansion (rather than proceeding iteratively, order by order), 
and to the definition of the resummed variables d and d that turn out to coincide with the standard first-order 
potentials in Poisson gauge. 

We now briefly consider the linear vector sector. Again defining a resummed variable, let us consider the vectorial 
potential uif. 


<+ = B f + -Bf 


(7.7) 


Linearizing Eqs. (14.131) and (14.91) . we immediately obtain 


^7 VV = ~^^Ga 2 p (V - , 

(7.8) 

2^(V • * (V • - 0. 

(7.9) 


Therefore, the linearization of the 1PF approximation also reproduces the equations of linear relativistic perturbations 
theory in the vector sector. Note that Eq. (17.81) is the first-order analogue of the nonlinear Newtonian Eq. (I6.3bl) . Eq. 
m instead gives the evolution of the first-order metric vector potential, which decay. In the nonlinear case, it is 
indeed shell-crossing that generates vorticity [90]. In general, nonlinearity generates the vector modes in the energy 
current that sources the gravitomagnetic frame-dragging potential jl, 2j. 

Finally, a remark on the tensor sector. Had we extended our post-Friedmann scheme to the 2PF order, in the linear 
regime we would have recovered the wave equation for free gravitational waves. In the nonlinear regime gravitational 
radiation would be generated by nonlinear sources terms, of higher order in scalar and vector modes, much in the 
same way that is generated at second order in standard perturbation theory []j} HJ [96j, [97} ■ 


VIII. NONLINEAR SET OF 1PF EQUATIONS 
A. Resummed metric potentials and new conservation equations 


We now recast the nonlinear equations obtained in the previous sections, using a suitable change of variables. In 
the previous section we have introduced the potentials (pp and ipp: we have shown that they coincide in both the 
linear and the Newtonian regime. It is then convenient to consider the following combinations: 

<Pg ■= -^{(Pp + 4>p), (8.1) 

^2 D p := 2 ( 8 - 2 ) 

the first generalizes the definition of the gravitational potential given in Sec. ED equation (16.81) . the second defines a 
new nonlinear post-Friedmannian quantity, Dp, negligible in the Newtonian and in the linear regimes 4 . With these 
new variables, up to 0(l/c 4 ), the metric can be written as 


ds 2 = —c 




dt 2 — 2 a-f-dtdx 1 + a 2 
c z 


<Pg . <Pg , D p 


1 + 2 ——fr + 2-y- 


Sij H —jhij ) dx l dx° 


(8.3) 


Moreover, it is convenient to define a new velocity variable 


* 1 
v i = v i -> 


(8.4) 


A difference in the two scalar potentials is expected from nonlinearity in GR, consistently with second-order perturbation theory results 

oim& 


4 
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representing the velocity of matter with respect to observers moving along the normal to the slicing. With these new 
variables, Eqs. (TOT) . (TOl) . m and 63} become 


g ° o + a = ^t° 0 


G\ + A6j = ^T\ 


G'. + Ad' = 


r-4 ‘ 




- tracefree 


4v 2 ^ - 4 


V 2 D p + 3a 2 ( H ) - V 2 4 + 


= -^t4tt Ga 2 pS + ^jAnGa 2 p(l + 8)v* 2 , 


4 S7 2 Dp + (j)G,k(t>G k + 6° 2 ( 4 -</>G + 2 ~(f>G + ( - ) </>G + </>G 

\ a a \a' 


1 


= —r 87 rGa 2 p(l + S)v* 2 , 


4 ( 2 ( D p,’i ~ + 2 0G,i<^ G J - \<t>G,k(t>G k S i 


(~V' • ^ (V • 


:V 2 # 


= -487rGa 2 p(l + <5) (4>*^ - ^u * 2 6j ) , 


(8.5a) 


(8.5b) 


(8.5c) 


o _ 8vrG 0 

^ » „4 1 i 


— V 2 Wi + 2 — (f>G,i + 2 </>G,i 
2 a a 


-r ■! — 2—Dp i + 2,D Pi 


1 

2 a . 


4a (j)G4>G,i + 2<l>G<j>G,i) + 8d(j)G4>G t i + 2u >k,i4>G,k ~ WjV 2 0g ~ 2Wfc0G,fei 

(8.5d) 


= —^ 87 rGap(l + S)v* - ^8nGap(l + 8) \y* (v* 2 — 4</>g)] ■ 


The continuity equation (15.71) and the Euler equation (15.61) become 


dS v*\, r , 1 


' dt a ala 


= 0 . 


dv* d * 1 1 

~T~ + ~V i + ~(f>G,i + — 

dt. a a c z 


-<t>G i(4^G + 'a* 2 ) ^ Stf-TT + —Dpi - -v*Vj(t> G ’ 3 - -V* 2 V* + -Uj,iV* J + -u 3 v*- 
a dt a a J a a a 


( 8 . 6 ) 


= 0 


(8.7) 


where here and in the following the convective derivative of a quantity Q is defined in terms of v*: 


dQ_ _ a . vlQ,i 

dt a 


( 8 . 8 ) 


B. Constraint type and evolution type equations 


Starting from the new definitions given above, we want to obtain a final set of equations, providing a unified 
description of gravitational instability in both the linear and the nonlinear regimes, valid for fluctuations on all scales. 

From Eas. (|8.5al) and ()8.5bl) we obtain an equation involving only the potential </>g, sourced by the matter and 
velocity perturbations: 


1 2 


rV 2 </>c 


4>g 


2 ~4>g 
a 


2 ~4>g — 
a 


4>g 


2 —,2 ±2 ^ 


3 V ^G - \4>G,i4>d l 


1 8 ttG 


a 2 pS 


+ —r4:TrGa 2 p(l + 8)v 


*2 


(8.9) 


Note that the linearized version of this equation is a combination of Eqs. (17. Gal) and (I7.6bl) . 
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A constraint equation for Dp is obtained from Eq. (14.121) [or, equivalently, by applying the operator djd 1 on both 
sides of Eq. (18.5cl) ]: 


1 2 


V 2 V 2 D/ 





,1 


.3 




J- 47 T Ga 2 p 


(1 + <J) 


„*2 


si 


( 8 . 10 ) 


Rewriting the definitions of Aj and Sj, Eqs. (14.161) and (14.171) . in terms of 4>g and v* 1 , also neglecting a higher-order 
contribution for Dp, we obtain 


Al = 2 <f>G,i<j>G ~ ^i4>G,k<t> G , 

S{ = (1 + 8){v*v* j - ^v* 2 Sj) . 


With these definitions Eq. (18.101) becomes 


11 , 

c 4 3 


2 V 2 Dp = - is nGaipS^ . 


k,j 


Finally, let us rewrite Eq. (14.181) as an evolution equation for the frame-dragging potential cj,; 


V 2 V 2 j — Wi H- U)i I — —T 


1 a 


2a 


a a 


h ( V2 A 


- 

J k,ji 


+ 


8ttGp (v 2 SI 


J °k,ji 


while the constraint Eq. (14.191) for hij remains formally unchanged, 


1 V 2 v 2 v 2 / l i = L 


-A l $ - V 2 A$6j + 2V 2 A-j + 2W 2 Al ki 6 lj - 2V 2 V 2 ^ 


+87rGa 2 p (-5® - + 2V 2 l S i fe fc ^ + 2V 2 5^ - 2V 2 V 2 Sf 


( 8 . 11 ) 

( 8 . 12 ) 


(8.13) 


(8.14) 


(8.15) 


At 1PF order the variables hij and Dp are not dynamical, i.e. they do not satisfy an evolution equation, rather they 
are given in terms of other variables by constraint equations. We can, therefore, conclude that the 1PF correction to 
Newtonian gravity introduces three nondynamical geometrical degrees of freedom (two d.o.f. in hij and 1 d.o.f. in 
Dp) and provides the evolution equation for the gravitational potential (f>G and the frame dragging vector potential 
Wj. For comparison, at OPF Newtonian order we only have three nondynamical d.o.f.: one in the scalar gravitational 
potential and two in the gravitomagnetic vector potential. In full relativistic theory we would have six dynamical 
d.o.f., as is already apparent in first-order perturbation theory. Therefore, in order to recover a fully dynamical theory 
we should extend our post-Friedmann scheme to the 2PF order, cf. [78]. However, the equations above, given that 
they extend the fully nonlinear OPF Newtonian equations and (as shown in Sec m include the scalar and vector 
first-order perturbation equations, are sufficient to study structure formation at all scales. 


C. New matter variables and simplified conservation equations 


There exists a further change of variable in the matter density, such that the 1PF continuity equation formally 
takes the usual Newtonian form [67]. 

Defining 


P ■= a 3 p{-g) 1/2 u° = p 





(8.16) 


and 8 := (p — p)/p, equation (18.61) becomes 

§ + —(8+ 1) = 0 ; (8.17) 

therefore “the mass defined in terms of the density p is conserved” |67] in the usual Newtonian sense 5 . Similarly, we 
can introduce a new velocity field v-i , such that the 1PF Euler equation simplifies. Defining 


1 + ( )!:v* 2 - 3</> G 


(8.18) 


5 Of course, mass is covariantly conserved in GR SB- For a discussion in second-order perturbation theory see [9^1 . 
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equation (18.71) takes the form 


dvj 

dt 


a „ 1 1 

H Vi H 0G,j "I - n 

a a c z 


4>G,i 

a 


i^+r 2 ) 


Dpi 

a 


+ 


Uj,jV J 


= 0 . 


(8.19) 


Let us note, however, that using this definition of velocity in the continuity equation (18.171) we would get some extra 
terms. 


IX. CONCLUSIONS 

Nonlinear structure formation in cosmology is traditionally studied with Newtonian N-body simulations and various 
approximation methods, e.g. Lagrangian perturbation theory. Relativistic perturbation theory and other approxima¬ 
tions such as a gradient expansion are used to study small fluctuations in the early Universe, in the CMB and on very 
large scales in the matter era. Thus there is a gap between methods used to study large and small scales. Various 
authors fS. 2t Ti-i60l| have recently pointed out that, given the high precision of current and future galaxy surveys, it 
is timely to investigate possible GR effects on structure formation at all scales. 

Assuming a standard flat ACDM cosmology, and a fluid approximation, in this Paper (Sec. IVUID we have developed 
a unified resummed nonlinear post-Friedmann formalism to study structure formation in the Universe. This relativistic 
scheme reduces to fully nonlinear Newtonian cosmology at leading order (Sec. IVII) and, if linearized (Sec. IVIII) . to 
standard first-order relativistic perturbation theory. Thus our post-Friedmann formalism is valid on all scales, bridging 
the existing gap in current studies of nonlinear structure formation. 

We have focused on obtaining a set of approximate nonlinear equation, consistently using the complete set of Einstein 
equations and the conservation equations. Rather than using the standard iterative procedure of the perturbative 
analysis where the equations are derived and solved order by order, we have constructed a resummed approximation 
scheme, which includes the first relativistic corrections to the equations of Newtonian cosmology, where a set of 
appropriately defined resummed variables satisfies a set of nonlinear equations. 

Our equations could be used to implement GR corrections in N-body simulations. This would be important 
in order to take into account causal, retardation and other GR effects that may be non-negligible for simulations 
aiming at a 1% accuracy [H] on scales of the order of the Hubble horizon. Indeed, we should consider that the 
relevant fluctuations for the formation of large-scale structure have not always been much smaller than the Hubble 
scale in the past. For instance, the present physical length of the horizon scale at decoupling is of the order of 
cf( 2 ec (l + Zdec) ~ 80ft -1 Mpc. Since galaxy surveys are going to cover a large proportion of the Hubble volume and aim 
at high precision measurements [42|, we should consider GR corrections for the large and intermediate scales where 
the Newtonian approximation is not good enough. A first step in this direction has been done in jJ. [2||. computing 
the frame-dragging gravitomagnetic vector potentials from N-body Newtonian simulations, and in [88|, where weak- 
lensing has been considered in the post-Friedmann framework. Our equations could also be used to include GR 
corrections in approximate Newtonian studies of various nonlinear effects such as BAO and CDM halos fl6l - [2ll |. A 
numerical implementation of our resummed equations is now a timely as well as doable goal: codes going beyond 
the quasistatic approximation have now been successfully developed for modified gravity [99j |. and codes specifically 
including GR corrections are in the making (56j . In this paper we have provided the theoretical framework aimed at 
this goal. 
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Appendix A: Some Useful Quantities 

In this Appendix, for completeness we collect some useful expression for the energy-momentum tensor and some 
geometrical quantities. 


Too = pc 2 + P {v 1 2 3 - 2U N ) + -^[v 4 5 6 7 8 ~ 4U P + 2 v 2 V n + 2 U N 2 ] 

T 0i = -paevi + — [B N i - Vi(v 2 + 2V n )\ 
c 


Tij = pa ViVj + 


2 

J+ C 2 


(4 Vn -b 2 L^v + v )viVj 2 B uVj) 


(Ala) 

(Alb) 

(Ale) 


Riemann Tensor: 


Roo :i 2 \ccSjj T Upf,ij\ T 
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R< 


+ aB N y^ + dB N ^j^ + 2Up t ij — Un^Unj — 2Um,(Yn,j) — 2?7jvC^jv,jj] 
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Oyfc - c 3 


2a ( V/vqfc -I f/jv,[fc I <5j]i + aB y. 
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Ricci Tensor: 
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Ricci scalar: 
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